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GLOBAL CLASSICAL SOLUTIONS OF THE “ONE AND ONE-HALF” 
DIMENSIONAL VLASOV-MAXWELL-FOKKER-PLANCK SYSTEM* 

STEPHEN PANKAVICH t AND JACK SCHAEFFER* 


Abstract. We study the “one and one-half’ dimensional Vlasov-Maxwell-Fokker-Planck system 
and obtain the first results concerning well-posedness of solutions. Specifically, we prove the global- 


in-time existence and uniqueness in the large of classical solutions to the Cauchy problem and a gain 


in regularity of the distribution function in its momentum argument. 
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1. Introduction 

From a mathematical perspective, the fundamental non-relativistic equations 
which describe the time evolution of a collisionless plasma are given by the Vlasov- 
Maxwell system: 


dtf + v-y,f + {E + vxB)- V„/ = 0 
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(VM) 


dtE = W X B — j, V ■ E = p 
dtB = -V xE, V • R = 0. 


Here, / represents the distribution of (positively-charged) ions in the plasma, while 
p and j are the charge and current density, and E and B represent electric and 
magnetic fields generated by the charge and current. The independent variables, 
t > 0 and x,v € represent time, position, and momentum, respectively, and 
physical constants, such as the charge and mass of particles, as well as the speed of 
light, have been normalized to one. 

In order to include collisions of particles with a background medium in the physical 
formulation, a diffusive term is added to the Vlasov equation in (VM). With this, the 
equations are referred to as the Vlasov-Maxwell-Fokker-Planck system. Since basic 
questions of well-posedness remain unknown even in lower dimensions, we study a 
dimensionally-reduced version of this model for which a; £ M and £ R^, the so- 
called “one and one-half dimensional” analogue, given by 


dtf + Vldxf -b K ■ Vyf = Ayf 

Ki = El + V2B, K2 = E2 — viB 
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(VMFP) 


dtE2 = -dxB - j2, dtB = -dxE2, d^Ei = p, dtEi = -ji. 


This system is the lowest-dimensional analogue that one may study and include elec¬ 
tromagnetic effects, as imposing r; £ R changes the model into the one-dimensional 
Vlasov-Poisson system. In (VMFP) we assume a single species of particles described 
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by f{t, X, v) in the presence of a given, fixed background (f) € (R) fl (R) n (R) 

that is neutralizing in the sense that 


4>{x) dx 


/(O, X, v) dv dx. 


The electric and magnetic fields are given by x) = {Ei{t,x), E 2 {t,x)) a.nd B{t,x), 
respectively. For initial data we take a nonnegative particle density with bounded 
moments G L^(R^), along with fields E^jB^ G Additionally, we 

specify data for Ei, namely 

{EiDAT) Ei{0,x)=J f{y,w)dw - 4>{y'^ dy. 

In fact, this particular choice of data for Ei is the only one which leads to a solution 
possessing finite energy (see [5] and [E]). The inclusion of the neutralizing density </) 
is also necessary in order to arrive at finite energy solutions for (VMFP) with a single 
species of ion. 

The analysis of (VM) has seen some progress in recent decades. For instance, the 
global existence of weak solutions, which also holds for the relativistic system (RVM), 
was shown in [3]. Unlike its relativistic analogue, however, no results currently exist 
that ensure global existence of classical solutions. Hence, the current work is focused 
in this direction. Alternatively, a wide array of results have been obtained for the 
electrostatic simplification of (VM) ~ the Vlasov-Poisson system, obtained by taking 
B = 0 within the model. The Vlasov-Poisson system does not include magnetic 
effects, and the electric field is given by an elliptic equation rather than a system of 
hyperbolic PDFs. This simplification has led to a great deal of progress concerning 
the electrostatic system, including theorems regarding the well-posedness of solutions 
[lOl im [141 US] • The book can provide a general reference to information concerning 
kinetic equations of plasma dynamics, including (VM) and (VMFP). 

Independent of these advances, many of the most basic existence and regularity 
questions remain unsolved for (VMFP). For much of the existence theory for collision¬ 
less models, one is mainly focused on bounding the velocity support of the distribution 
function /, assuming that /° possesses compact momentum support, as this condition 
has been shown to imply global existence [7]. Hence, one of the main difficulties which 
arises for (VMFP) is the introduction of particles that are propagated with arbitrarily 
large momenta, stemming from the inclusion of the diffusive Fokker-Planck operator. 
Thus, the momentum support is necessarily unbounded and many known tools are 
unavailable. Though the x-support of the distribution function is not bounded, we are 
able to overcome this issue by controlling large enough moments of the distribution 
to guarantee sufficient decay of / in its momentum argument. This also allows us 
to control nonlinear terms that arise within derivative estimates. As an additional 
difference arising from the Fokker-Planck operator, we note that when studying colli¬ 
sionless systems, in which Ayf is omitted, L°° is typically the proper space in which 
to estimate both the particle distribution and the fields. With the addition of the 
diffusion operator, though, the natural space in which to estimate / is now L^. Thus, 
to take advantage of the gain in regularity that should result from the Fokker-Planck 
term, we iterate in a weighted setting. Other crucial features which appear in¬ 
clude conservation of mass, and the symmetry of the diffusive operator. The main 
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advantage of the diffusion operator is that it allows one to estimate spatial deriva¬ 
tives of the density in independent of the momentum derivatives. This is not 

true for the Vlasov-Maxwell system, which is conservative rather than dissipative. 
Additionally, the appearance of the Laplacian allows the particle distribution to gain 
regularity in its momentum argument in comparison to its initial data. Finally, we 
note that our methods utilize an extra conservation law arising from the structure 
of the one-and-one-half dimensional system in order to bound the electric and mag¬ 
netic fields. Hence, they do not immediately apply to higher-dimensional analogues 
of (VMFP), though many of the other ideas presented below will likely be useful in 
the two, two-and-one-half, and three dimensional settings. 

Though this is the first investigation of the well-posedness of (VMFP) in the large, 
others have studied Vlasov-Maxwell models incorporating a Fokker-Planck term for 
small initial data. Both Yu and Yang m and Chae [1] constructed global classical 
solutions to the three-dimensional Vlasov-Maxwell-Fokker-Planck system for initial 
data sufficiently close to Maxwellian using Kawashima estimates and the well-known 
energy method. Additionally, Lai arrived at a similar result for a one and one- 
half dimensional “relativistic” Vlasov-Maxwell-Fokker-Planck system using classical 
estimates. The system in this work features a relativistic transport term, but still 
utilizes the Laplacian A„ as the Fokker-Planck term. We note that the relativistic 
transport operator yields an extremely beneficial result, known as the cone estimate 
(see 0), whereas the non-relativistic transport within (VMFP) does not. Thus, one 
essential component of the current paper is to overcome the lack of bounds on energy 
inside the light cone. Finally, we mention |12j . which arrived at similar results to 
our own but studied the relativistic Vlasov-Maxwell system with a Lorentz-invariant 
diffusion operator. While we utilize some of the tools introduced within and 
related articles 01 [la, we also introduce a number of new methods to overcome the 
loss of the cone estimate, finite speed of propagation, and a priori field bounds in 
order to arrive at the first large data global classical solutions to (VMFP) set in any 
dimension, see Theorem 1.2 below. First we state a local existence theorem: 
Theorem 1.1. Let a > 8 and e > 0 and denote 


Vo = + kP- 

Assume that cj) G C'^(R)niL^(M)nL^(M). Assume that is continuous, nonnegative, 
and bounded and possesses a partial derivative with respect to x such that 


jj dv dx + jj Up ^+‘"(5x/°)^dudx 

is finite. Assume that £’ 2 , G C^(M) fl Then there zs T > 0 depending only 

on 


II dvdx+\\E°\\j,^ + \\BYm, 

f G C'([0,r] X R^) nC'^((0,T] x R^) with second order partial derivatives with respect 
to vi,V 2 that are continuous on (0,T] x R^, and {E,B) G C'^([0,r] x R) for which 
(VMFP) holds, [EiDAT) holds, and 
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Moreover, f is nonnegative and bounded, and 


I[ dvdx+\\Emm + rnmn^ 

is bounded on [0,T]. Lastly, the above solution is unique. 

Note that /° is not assumed to be smooth in v. Now we may state the main 
result: 

Theorem 1.2. In addition to the hypotheses of Theorem 1.1, assume that E 2 ,B^ € 
L^(]R) and Vq/^ G for some S > a + 2 + s, and v^f^ G Then, the 

local solution of Theorem 1.1 may be extended to [0, c») x R^. 

We note that Theorems 1.1 and 1.2 can be altered to accommodate a friction 
term. In the model with friction, the Vlasov equation is changed to 


dtf + vA + fc • V,/ = V, • (V„/ + vf). 

The new term is lower order and does not change either of the results. 

As additional evidence of the gain in regularity in v we also state: 
Proposition 1.3. Assume the hypotheses of Theorem 1.2 hold. Then for all t > 0 

JJ + 1 |V„/|' + |V2/|"^ dvdx < Ct. 


This paper proceeds as follows. The proof of Theorem 1.1 is postponed to Section 
4 and Sections 2 and 3 assume the result of this theorem. In Section 2 we state six 
lemmas and show how Theorem 1.2 follows from them. The proofs of these lemmas 
and Proposition 1.3 are contained within Section 3. 

Throughout the paper C denotes a positive generic constant that may change 
from line to line. When necessary, we will specifically identify the quantities upon 
which C may depend. Regarding norms, we will abuse notation and allow the reader 
to differentiate certain norms via context. For instance ||/(t)||oo = sup \f{t,x,v)\ 

whereas ||i?(t)||oo = sup x)|, with analogous statements for || • ||2 and < •,• > 

which denote the norm and inner product, respectively. 

2. Global Existence 

Throughout this section we assume the hypotheses of Theorem 1.1 hold. Let T be 
the maximal time of existence and, in order to prove Theorem 1.2 by contradiction, 
assume T is finite. 

To begin, we will first prove a result that will allow us to estimate the particle den¬ 
sity and its moments. When studying collisionless kinetic equations, one often wishes 
to integrate along the Vlasov characteristics in order to derive estimates. However, 
the appearance of the Fokker-Planck term changes the structure of the operator in 
(VMFP), and the values of the distribution function are not conserved along such 
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curves. Hence, the following lemma (similar to that of [2]) will be utilized to estimate 
the particle distribution in such situations. 

Lemma 2.1. Let g G L^((0,T), and ho G n be given. Let 

F{t, X, v) = F{t, X, v) + B{t, x){v 2 , —vi) be given with F G T) x R^; R^) and 

B G W^’°°((0,T) X R;R). Assume h(t,x,v) is a weak solution of 


Lh 


dth + vidxh + F{t, X, v) • Vyh — Ayh = g{t, x, v) 
h(0, X, v) = h^ix, v) 


so that h G L^((0,r) X satisfies 


( 2 . 1 ) 



h {-dtfi - vidxfi) + ^yh ■ {F(j) + Vyfi) - gfi 

~JJ ho{x,v)(j){0,x,v)dvdx = 0 


dvdxdt 


for every (j) G 2?([0,T) x R^). Then, for every t G [0,T] 


||/l(OI|oo < ll/lolloo + / ||g(s)||oods. 

Jo 

Another useful tool will be the conservation of mass and energy growth identities, 
which we establish in the next result. 

Lemma 2.2. (Conservation Laws). Assume v^f^ G L^(R^). Then, for every t G 
10, T), 


ll/Wlli = ll/°lli 

and 


\v\‘^ f{t,x,v) dv dx + / + B‘^)dx < C{1 +t). 


Next, we state a lemma that will allow us to control moments of the particle 
distribution. 

Lemma 2.3. (Propagation of t! 2 -moments). Let p G [0, oo) be given and assume the 
hypotheses of Lemma 2.2 with E 2 , G L^(R). Let R{s) = Vl + . If\\R{v 2 )^f^\\oo <| 

00 , then for any t G [0, T) 


\\R{v2rfmo. < ct. 


With control of velocities in the V 2 direction, we are able to control the induced 
electric and magnetic fields. Bounds on moments of the particle density then follow 
from this result. 

Lemma 2.4. (Control of fields and moments). Assume there is 6 > 4, such that 
v^f^ G L°“(R^),Uq/° G L^(R^), and B^ G L^(R). Then, for any t G [0,T) 
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lkVWI|oo<CT 


( 2 . 2 ) 


\\E{t)\\^ + \\B{t)\\^<CT 


(2.3) 


and 


Vq ^f{t) dv < Ct 


(2.4) 


OO 


for any /3 G [0, ^). 

Thus, once control of the fields is obtained, any higher moment of the particle 
distribution function can be controlled as well, assuming that the initial distribution 
possesses the same property. Next, we utilize energy estimates to bound the density 
and its derivatives in L^(K.^). 

Lemma 2.5. Assume the hypotheses of Lemma 2-4 hold, then for every t G (0,r] 


j^\\fmi = -n^vfit)\\i 


and thus 


ll/WI| 2 <||/°|| 2 . 

If additionally, G L^(K^) for some 7 > 0, then 



and thus 


\Hfm2 < Ct 


for every t G [0, T). 

Lemma 2.6. Assume the hypotheses of Lemma 2.4 hold with 6 > 8. Then for every 



lk?5,/(t)|U2 + \\d,,E{t)\\L^ + \\d^B{t)\\L^ < Ct. 


Proof. Now we may prove Theorem 1.2. Applying Lemma 2.5 with 7 


2 


yields 


// 



dv dx < Ct- 
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Applying Lemma 2.6 with 7 


CL — 2 -j- £■ 
2 


yields 


JJ v^-^+^id,ffdvdx + Ji\d,Ef + {d,Bf)dx<CT. 

Also by Lemma 2.2 


J {\E\^ + B^)dx<CT. 

Taking (/(t), E 2 {t), B{t)) as an initial condition and applying Theorem 1.1 we find the 
solution may be extended to [ 0 , t + r] with r > Ct- This contradicts the maximality 
of T and completes the proof. □ 

3. Proofs of Lemmas and Estimates The first result (Lemma 2.1) is very 
close to a previous lemma m, in which this property was shown for the relativistic 
Fokker-Planck operator. One alteration necessary in the proof of [12, Lemma 1] is to 
change the relativistic velocity vi to vi, which does not affect the conclusion. Also, 
here E is not in L°°, but Wy ■ F = Wy ■ F and the proof of [12, Lemma 1] still applies. 
Hence, we omit any additional details. 

Proof. [Lemma 2.2] We begin with conservation of mass. Integrating the Vlasov 
equation over all (x, v) we find 


d 

dt 



f{t, X, v) dv dx 


= 0 . 


Thus, using the decay of /° we find for every t £ [0,T) 


Jj 


') dvdx 


/°(x, v) dv dx < 00 . 


(3.1) 


To arrive at the estimate of the total energy, we multiply the Vlasov equation by 
|up and integrate in v. The Fokker-Planck term becomes 


J \v\^Ayfdv = —J 2 u-V„/du = 4 J f 

after two integrations by parts. Hence, using the divergence structure of the Vlasov 
equation, we arrive at the local energy identity 


where 


dte + dxTn = 4 


f{t,x,v)dv 


(3.2) 


e{t,x) = J \v\^f{t,x,v) dv + {\E{t,x)\'^ + \B{t,x)\'^) 


and 
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m{t,x) = J Vi\v\'^f{t,x,v)dv + 2E2{t,x)B{t,x). 
We integrate (13.21) over all space to deduce the global energy identity 


d 

dt 


J e{t,x)dx = 4: jj f^{x,v)dxdv 


whence we find 


J e{t, x)dx < C{1 + t) 

for all t G [0, T). □ 

Now we utilize the conservation laws to prove Lemma 2.3. 

Proof. [Lemma 2.3] We begin by bounding the potential associated to the electric 
and magnetic fields. By Lemma 2.2 we have 


J \j 2 it,x)\dx < C{l+t) 


and hence 


J J j2(T,y±{t-T))dTdy < J \j2{T,y ± {t - t))\ dydr < C{1 + tf 
Since G it follows that 


J \B{t, a;)| dx < C{1 + 


and we may define 


A{t,x) = f B{t,y)dy. 

J —OO 

Note that dxA = B and dtA = —£ 2 - Moreover, using Maxwell’s equations, we find 


{dt-dt)A = j2 


and thus 


A{t,x) =-{A{0,x - t) + A{0,x+ t)) + - / V2f{s,y,v)dvdyds. (3.3) 

t/0 'J X —t-j-s J 


The (x, u)-integral can be bounded using Cauchy-Schwarz and Lemma 2.2 as 




J j V 2 f{s,y,v)dvdy < f{x,y,v)dvdy^ vlf{x,y,v)dvdy 

< \\f\\l^^(^JJ\v\^f{s,y,v)dvdy 


1/2 


\ 1/2 


< cii + sY/"^. 

Hence, using the assumptions on initial data and integrating, we find 

PWlloo <C(l + t)3/2 <Cr. 

Next, we utilize the identity 


dtA + vidxA = —E2 + viB = —K2 

within the Vlasov-Fokker-Planck equation. In particular, let tp £ C'^(K) be given and 
multiply this equation by 'ip{v 2 + A{t,x)). Denoting the VFP operator by 

Vh := dth + vidxh + K ■ Vyh — Ayh, 

we find 


V(V'(v2 + A)f) = -2'ip'{v2 + A)dy^f - 'ip"{v2 + A)f. (3.4) 

Next, define the function R{x) = vTT^. To prove the first assertion, we use 
'ip{x) = R’^{x) within (13.41) and derive the equation 

V{RP{v2 + A)f) = -2p{v2 + A)RP~'^{v2+A)dy^f-pRP~'^{v2+A)[l + {p-l)\v2 + A\'^]f. 
Using the identity 


dy^f = R P{v2 +A)dy^{RP{v2 +A)f)-p{v2+A)R '^{v2 + A)f 


the right side becomes 


-2p{v2+ A)R '^{v2 + A)dy^{RP{v2+A)f)+pRP {v2 + A)[-1 + {p + l)\v2 + A\'^]f. 
Hence, if this first term is included within the VFP operator by defining 


K = K + 



V2 + A \ 

1 + |i ;2 + Hp / 


to form the new operator V, we find 
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ViRP{v2+A)f)=pRP-\v2+A)[-l + {p + l)\v2+A\'^]f. 

We note that the term on the right side satisfies 


\pRP-Hv 2 + A)[-l + {p+ l)|r ;2 + A\^]f\ < CRP-\v 2 + A)f < CRP{v 2 + A)f. 
We invoke Lemma 2.1 with h = RP{v2 + A) f and £ = V so that 

\\RP{v 2 + A{t))f{t)\\oo < \\RP{v2 + ^( 0 ))/°||oo + c f \\CRPiv2 + A(s))/(s)||oo ds. 

Jo 

By Gronwall’s inequality we find 


\\RP{v2+A{t))f{t)\\^<CT 

for t G [0,T). Finally, the previously established control of ||^(t)||oo yields the first 
result as for p > 0 


R^{v2)fit,x,v) = {1 + \V2 + A{t,x) - A{t,x)\'^)P^'^f{t,x,v) 


< C{Rp{v 2 + A{t,x)) + \A{t,x)\P)f{t,x,v) 


< c\mv2+A{t))f{t)\\^ + \\Amafm 


< Ct- 

Hence, taking supremums we find 


\\R^{v2)f{t)\\^<CT. 


□ 

Using this result, we may bound the fields and moments of the distribution func¬ 
tion. 

Proof. [Lemma 2.4] We first bound Ei using conservation of mass so that 


ll-E'i(i)||oo = sup 
xeM. 


J (^J f{t,x,v)dv — (j){x)^ dx < JJ f{t, x,v) dvdx + {{(pWl < C 


Next, we estimate the other field components. Using the transported field equations, 
we find 


{E2 ± B){t, x) = {E2 ± B){ 0 ,x ^ t) — J J V2f{s, x =F (t — s), v) dvds. (3.5) 

Note that E^jB^ G L°°{M.) by the Sobolev embedding theorem. Thus, for any ei, £2 > 
0 we have 
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\{E 2 ± B){t,x)\ < C{^l + J Jr {s,x T {t - s),v) 

\v2\R~'^‘^'^^'^\v2) ^(S) T (i — s), f;) dvd^ 

/ 1 + g'l 7— (l + gj) \ 

< C(l +j^\\R^{v^)f{s)U^ \\R‘i{v2)f{s)U ^ dsj 

(2 + £ 2)7 

where q = - - -r. We choose 7 > 1 + ei and note that <5 > 4 ensures that we 

7 — (1 + £ 1 ) 

may also choose q < ^ < S. Define the function 


F{t) := sup ||u;]'/(s)||oo 

sG[0,t] 

Invoking Lemma 2.3 with p = q we find 


\\iE2±B){t)\\^ < 


Ct 



sup ||i?'>'(-i;i)/(s 

sG[0,i] 



/ 1+^1 \ 
< CT[lEF{t) — y 


Using the identity 


(3.6) 


E 2 {t,x) = ^{[E 2 {t,x) + B{t,x)] + [E 2 it,x) - B{t,x)]) 

we see that the same bound holds for ||U 2 (t)||oo- 

Next, we multiply VFP by Vq and use the same method as in the proof of Lemma 
2.3 to derive the equation 


= ^iyE)f-2jVQ^v-Wy{v^f)+j{-j\v\'^-2)v^ ^f. (3.7) 

If the second term on the right side is included within the VFP operator by defining 


to form the new operator V, we find 

V(uJ/) = 'yv^~‘^{v E)f+ -f{j\v\^-2)v^~'^f 

=: I+ 11. 


Clearly, 
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II < c\\vr"m\\^ < cwv^ofmoo < cF{t). 

Estimating I requires the field estimates, which yield 

I < CuJ-"(|ui| + |z;2|-p2(t)||oo)/ 

< Ct (^F{t) + Fit)"^ + 

since 7/2 < S. We combine these estimates and invoke Lemma 2.1 with h = v^f and 
£ = V so that 


IkoVWIloo < |koV°lloo + Ct £ [Fit) + F{t)^ + F{s) 


'y-l + ei 


ds 


Taking the supremum in t and choosing ei < 1 


F{t) < F{0 )+Ct (e(s)+F( s)"^+E(s)^-^) ds < F{0 )+Ct {1+F{s)) ds 

Gronwall’s inequality then yields the bound F{t) < Ct for any t G [0, T) and 1 < 7 < 
6. The bound on moments of the distribution function follows immediately and the 
field bound 


||£2(t)||oo + P(t)||oo<GT 

then follows from (13.61) . Finally, using the bound on moments of the density, control 
of the u-integral follows since we have 


j 1^0 <\\vof{t)\\^- J ^^dv<CT 

for 7 < (5 and taking the supremum in x yields (1^ . 0 

Proof. [Lemma 2.5] We proceed by using dissipative estimates. First, we compute: 




{-vAf-K-y,f + A,fJ) 


= -(uia,/,/)-(K-V,/,/) + (A„/,/). 


Notice that the first two terms are pure derivatives in x and v, respectively. Thus, 


{vAf, f) = \ 


dx{vif^) dvdx = 0 
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and 


{K ■ V„/, /) = i yy v„ • {Kf) dvdx = 0. 

Finally (A^/,/) = -||V^/(t)|||. Hence 

and the first conclusion follows. 

Similarly, we may multiply by and proceed in the same manner. Within this 
estimate we will use uq > 1 in order to increase moments of the estimates where 
necessary so as to match the results of the lemma. Computing the time derivative 

^^Iko/Wlli = JJ f[-vid^f - K - Vyf + Ayf]dvdx 

= / + // + III. 

The first term vanishes as it is a pure ^-derivative. For II, we integrate by parts and 
use the field bounds of Lemma 2.4 so that 

II = -^JJvl'^V,-{Kf)dvdx 

= 7 JJ Vq'’~‘^v ■ Ef^ dvdx 
< CrHfml 

To estimate HI, we integrate by parts twice in the first term and once in the second 
term to find 


III = - JJv,{vly)-V,fdvdx 

= - JJ + VQ'^Vyf)-Wyf dvdx 

< C'lko/Wlli - lkdV„/(t)||i. 

Combining the estimates, we find 

< crKiml - 

as in the statement of the lemma. Additionally, because the second term on the right 
side of the inequality is nonpositive, we invoke Gronwall’s inequality and find 
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Kfml<CTKf\\l<CT 


for every 7 > 0 for which the norm of the initial data Ur'o/^lh is finite. □ 
Proof. [Lemma 2.6] If / were we could compute 


d 

dt 


^ {d^ff dvdx = 


jj - vidlf 


-K ■ Vydxf - dxK ■ Vvf) dvdx 

-2 JJ vl''^\Wyd^f\'^ dvdx + JJ {d:cf)'^{Ayvl'^ + K ■Wyvl'^)dvdx 


+2 / / dvdx 


27 \ 


< “2 yy VQ^\Vvdxf\ dvdx+ C JJ (dxf) VQ'^{l + \E\)dvdx 

+cj JJ f^\dxK\‘^VQ'^ dvdx JJ v^'^lVyOxfl'^ dvdx 

+ J JJ vl'^{dxff dvdx'^ . 


Using the inequalities —x^ + Ax < —A^ and 2xy < x'^ + yields 


^ / / vl"^ {dxfY dvdx < C 11 {dxf)‘^vl'^ {1 + \E\) dvdx 


+C 11 PIOxEI^Vq'^ dvdx 


and hence 



{dxff dvdx 


< C + C 


+ |£'|) dvdxdr 


+C 



Vq'^ f^\dxK\^ dvdxdr. 


(3.8) 


By a standard regularization argument it follows that ()3.8I1 holds for the solution 
{f,E,B) with the regularity stated in Theorem 1.1. Applying (12.31) and (12.41) with 
/3 = 27 + 4 yields 
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JJ Vo'^idxf)^ dvdx < C + Ct J JJ vp{dxff dvdxdr 

+C m{\dxE\^ + {dx^Bf)\\f{t)\\L^ 

< C + Ct J JJ {dxf )^ dvdxdr 


dvdxdr 


+Ct I I i\dxE\^+ {dxBy)dxdT. 


Similarly, if E and B were we could compute 


(3.9) 


so 


^ J[\dxE\^ + {dxBfjdx = -2jdxE-dxjdx 

< 2^J\dxE\'^dx y Jldoojl'^dx 

J[\dxEf + {dxBf]dx < C + 2^ ^ J\dxE\^dx ^Jldx^jfdxdT 


(3.10) 


< I i\dxE\^ + \d,j\^)dxdT. 


By a standard regularization argument it follows that (13.101) holds for E and B € C^. 
Since 7 > 2 we have 

Idxjl"^ < (^J \dxf\vodv^ 

< J \dxf\^vl’^dv J vl~‘^"^dv<C J \dxf\^vl’^ dv 

so adding (13.91) and (13.101) yields 

JJ {dxf)"^ dvdx + J {\dxE\'^ + {dxB)'^) dx 

< C + Ct j (^J \dxE\'^ dxdr + JJ {dxf)"^ dvdx'^ dr. 

An application of Gronwall’s inequality completes the proof. □ 
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Proof. [Proposition 1.3] If / were we could compute the following: 


dt 

dt 




jj |V„/p dvdx 

II \Vlffdvdx 


= -2 


II |V^/|^dndx, 

-2 II + dyjda;f)dvdx, 


= -2 / / \Vlf\ dvdx - A // V vdv^f -V vdxf dvdx 


= -2 


II dvdx + AII dxfAydyJ dvdx, 


SO 




V^/pdwda:: —t // | V^/p dndx — 2t // dy^fdxf dvdx 


11 + AdxfAydyJ)dvdx 


< - II \yvf? dvdx + 2tJ (dxf)^ dvdx W {dy^f)'^dvdx 


+t \ - // dvdx+ 2 J (dxf) dvdx , {AydyJY dvdx 


Using the inequality — + Ax < -A^ twice yields 


^ //(f+t|V„/|2 + it2|v2/|2)dudx 


< / / {dxfY dvdx 


and 


(/'+i|V„/P + -t2|V2/|2)dudx 


< II(f^)^ dvdx + CI II{dxfy dvdxdr. 


(3.11) 


Again, by a standard regularization argument it follows that p.lll) holds for the 
solution constructed in Theorem 1.1. By Lemma 2.6 
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{dxff dvdx < Ct 


so the proposition follows from (13.111) . □ 

4. Local Existence 

Define 6 = a — 4, 


-F= /: 


[0, oo) 


v§f, v-^dxf 


and 


ll/lb = hSfWh-^ + hodxfWL-^. 

For the time being we consider smooth initial data (/°, i?°). 

Let T G (0,l),i? > I,-)/; : R —>■ [0,1] be smooth with s < —1 => ijj^s) = 1 and 
s > 0 => '(/'(■s) = 0. Define V'^(w) = ipi\v\ — R). For {E, B) G <^([0, T]; iL^(R)) smooth 
define 


by 


C^iE,B) = if,E,B) 


K = E + i/;^(n)B(v2, 


(4.1) 


dtf + v,dxf + K-Vj=Aj, /(0,.,.) = /O, (4.2) 


p = J fdv-(j), j = j vfdv, (4.3) 

El = f pdy, (4.4) 

J —OO 

dtE2 + dxB =-j2, dtB + dxE2=0, (4.5) 


iE2,B){0,-) = iE^„B°). (4.6) 

Note that K is bounded and hence by Proposition A.l of m, m has a solution 
/ G L^([0,T] X R;iJ^(R^)). Reference [H] may be used for this also. 

Let 0 = 0 + 2 +e, /3 = 5+ 2 + e, and 
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Co > \\EI\\h^ + WB^Whi + IHfh^ + 


We assume that 


|l(C,B)(t)||^i < lOCo (4.7) 

on [0,T]. Within the remainder of this section constants may depend on a,T and Cq 
but not on R or Vw/°. 

First, using (HTl) and the Sobolev embedding theorem, 


dt 


VQf^dvdx = — 2 / / VQlVyffdvdx 


+ f {AyVQ + E ■ V^Vq) dvdx 


(4.8) 


< C 11 dvdx. 


Hence, by Gronwall’s inequality 


VqP dvdx < C. 
Similarly, and using the Cauchy Schwartz inequality 

d 


dt 


vlo{dxf) dvdx = -2 vlQ\Vydj:f\^ dvdx 


+ 11 {dxf)^{/^vVQ + E ■ VyV^) dvdx 


+2 JJ fd^K ■ {v^V vdxf + d^fV yV^) dvdx 


< -2 


JJ v^\Vydxf\^ dvdx+ C JJ {dxffv^ dvdx 


(4.9) 


+C / / f{\dxE\ + \dxB\){v^+^\Wydxf\ + \dxM)dvdx 


< -2 11 v^lWyOxfl"^ dvdx + C hdxffv^ dvdx 


+C\I 11 pv^+\\dxE\'^ + {dxBY) dvdx 


Vol^vdxfl'^ dvdx 


v^{dxfY dvdx 
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(4.10) 








Using —x^ + Ax < and 2xy < x'^ + yields 


II i^oidxff dvdx < C 11 v^{dxff dvdx 


+C 11 pv^+^dv{\d,E\'^ + {d^B)'^) dx. 


(4.11) 


Also by (l4^ 


J Pvt^dv = lpVo^^'dv = 


^ + 0 


‘^fdxfvQ " dvdy 


< 2||/(i)uo" |lL2||a^/(i)i;2 11^2 


(4.12) 


< C + {dxfY dvdx 


so using (14.71) . (14.111) yields 


ilpSis-fy 


dvdx < C + C / / VQ^dxf)"^ dvdx. 


Hence 


jj vpdxff 


dvdx < C. 


(4.13) 


Next consider {£,B) S (^([O,T];i7^(R)) smooth for which (14.71) holds and define 
(F, f, B) = C^{£, B) where R <Tl and /C, F, P, J, £, and B are defined as in equations 
(14.11) through (I4.6|l . 

Let G = {E, B) — {£, B) and note that 


\K-IC\ < vo\G\+vo\B\\il;^{v)-il;'^{v)\ 

< vo\G\+v^+i\B\R-^ 


(4.14) 


and similarly 


\dxK-dxlC\<vo\dxG\+vp^\dxB\R-i. (4.15) 

Let g = f — P. Proceeding as before in (j4.10l) and (14.111) we have 
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11 Ka dvdx = 


-2 


JJ v^lVyg]^ dvdx 


+ g^{AyVQ + E ■ VyVQ) dvdx 


+2 / / F{K - K) ■ {v'^S/yg + gVyV^) dvdx 


< C JJ dvdx + C JJ F^\K - lC\^v^ dvdx. 

By (14.141) . the Sobolev embedding theorem, and (14.9|) we have 

F^\K - /Cp7;“ dvdx < [[ F^{G^ + B^R-^)v^ dvdx 


< ciwcmi^ + mmUR-n < + R-n- 

Substitution into (I4.16P and using Gronwall’s inequality yields 

JJ VQg^dvdx<C J |1 G(r)||^i dr + 

Again proceeding as in (14.101) and (14.111) we have 


_d 

dt 


J J Vo{do,gf dvdx =-2 JJ VQ\\7yd^g\dvdx 


+ (dxgf(AyVQ + K ■ Vyvlh) dvdx 


(4.16) 


(4.17) 


+2 / / {dxF{K - 1C) +gdxK + Fdx{K - K.)) ■ {vQ'C/ydxg + dxg'C/yVg) dvdx 


< C VQ^dxg)^ dvdx + C / / {dxF)^\K — JCfVgdvdx 


+C // g^\dxKfVQdvdx + C // F^\dxiK — ]C)fVQdvdx 


By (14.141) . (14.131) and the Sobolev embedding theorem we have 


{dxF)^\K - IC\^vj) dvdx 


< C {dxF)WGY + WR-^)vP dvdx 


(4.18) 


< G(|lG(t)||^, + \\B{t)\\]j,R--) < G(||G(t)||^, + R-^). 
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(4.19) 














Note that (using (14.171) ) 


< 2 


\g\\dxg\vQ^ dvdx 


< 2 


g^v^ dvdx 


(dxgfVQ dvdx 


< cj \\G{T)\\‘jjidT + CR + JJ{dxg)‘^VQdvdx 

so by (I4.20L (I4.7p . and the Sobolev embedding theorem 

f\dx:Kfv^odvdx< jj f{\dxE\^ + [dxBf)vl+‘^ dvdx 

< cJ ||G(r)||^i(ir + + C JJ{dxg)^VQdvdx. 

Using (I4.15L (14.121) . (I4.13p . and (14.71) we have 

F^\dx{K - K:)\^v^o dvdx 

< 11 F^{\dxG\^ + {dxBfR-^)v^ dvdx 


< G Ji\dxG\^ + {dxBfR-^)dx 


< C\\Git)\\l,+GR-^. 

Substitution of (14.191) . (I4.21L and (14.221) into (14.181) yields 


_d 

dt 


JJ VQ{dxg)^ dvdx < C JJ VQ^dxg)"^ dvdx 

+G||G(t)||^, + GR-^ + G f ||G(r)f^, dr. 

Jo 


By Gronwall’s inequality we have 


JJ VQ^dxgY dvdx < G J ||G(r)||^i dr + Gi? 


Next we consider the fields. We have 


— / {\E\^ + B‘^)dx = -2 I E- jdx. 
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(4.20) 


(4.21) 


(4.22) 


(4.23) 




















By (gH) 


/ 


\j\'^dx < 




dx < C 


so 


and 



dx < C 



_1 

2 


\E\^ + B^) dx < Cq + Ct 


(4.24) 


follows. Similarly by (14.131) 


dx < 


{dxffv^dv 



dx < C 


so 


J (^\dxEf + [dxBy^ dx < 2\\dxm\\L4dxjmL- 

< C\\dxE{t)\\L2 

and 


\dxE\‘^ + {dxBf] dx<Co + Ct 


(4.25) 


follows. In the same manner (14.171) yields 


\j - J\^dx < 




dx 


< 



\\G{T)\\jJ^dT + CR-H 


and, letting G = {E, B) — (f, S), 


\G\^dx<C \\G{t)\\]j^+GR-H 


follows. Lastly (14.231) yields 
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dx 


and 


\dx{j - J)\^ dx < 



{dxg) Vq dv 



< 



||G(r)||^i dr + CR-H 


dx<C ||G(r)||^idr + CR-H 


follows. 

From (I4.24|) and (14.251) we have 


\\iE,B){t)fH, <2Go + Gr< lOGo 
for T suitably restricted. Also, 


rt (4.26) 

< C \\iE,B)iT)-i£,B)iT)fH,dT + CR-H. 

Jo 

Define (/”+\ B”+i) = C^"{E'^,B'^) for n > 0 where f,E°,B° are deter¬ 

mined by the initial conditions. By (I4.26|) we have 

||(^n+l _ /n)(i)||2, + ||(F;"+l,i3"+l)(t) - {E^B^mrH, 

rt (4.27) 

< C ||(G",B")(r) - (F;"-\B"-1 )(t)||^i dr-f G2-’"^t. 

Jo 

Suppose A> 1 and 


0 < X^+\t) X^{T)dT + A-^t 


for n > 0. Then by induction 


‘(i) < 




nl 


{cty 


i\ 


< 




■A 


-n^CAT 


Since this bound is summable, it follows from (14.271) that (i?",i?") is Cauchy in 
G([0,T];i7i) and P is Cauchy in G([0,T];JF). Let if,E,B) = lim„^oo(r , A", B"). 
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We will now use the explicitly known fundamental solution for the linear equation 


dtf + vidxf = Ayf, 
namely, for 0 < t < t, a:, y G K, u, in G 


g{t,x,v,T,y,w) 


= [47r(t — r)] — 

- 3 


exp I —3 


t){vi + wi)f 


{t - rf 


This may be derived from line (2.5) of [16] by letting /3 —>■ O'*' (with N = 2) and then 
integrating in 1 / 2 - It may also be derived directly by Fourier transform. Since 

- K'- ■ V„/"+i 

,r+no,-,-) = /° 

it follows by Theorems II.2 and II.3 of [16] that 


=H+ II g{t, X, V, T, y, w){-K^ ■ V^r+i) 


dwdydr 


{'r,y,w) 


where 


H{t,x,v) = jj g{t,x,v,0,y,'w)f°{y,w)dwdy. 
It is easy to check that 


\'^wgit,x,v,T,y,w)\ dwdy < C{t - t) 2 


and it follows that 


/"+^ = H+ II W^g{t, X, V, T, y, w) ■ (iG"/"+') 


dwdydT. 


('r,V,w) 


By Lemma 2.1 


o<r<sup/° 


so 


0 </<sup/° 
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(4.28) 








follows. 

Thus far we have assumed {f'^, E 2 , B^) to be smooth. Now consider 
as in Theorem 1.1. Consider a sequence (/°^, 5°^) of smooth initial conditions 

with 


ll< (/°" - f)\\L^ + - /°)IIl= + ^ 0. 

By a limiting procedure like the above we conclude that (14.281) holds for (/°, 
as in Theorem 1.1. 

By Theorem 11.3 of [16] H G (^([OjOo) x R^) fl (^^((OjOo) x R^) with dvidv^H 
continuous on t > 0 and 


r dtH + vid,H = A,H 
\ = f. 

Next, we will show that / is differentiable in v and Holder continuous in x. 
Proceeding as in (4.16) we have 

Tt If < dvdx = -2 // ug I V„ (/"+! - /'+!) I" dvdx 

+ II (/”+' - ■ V,ug) dvdx 

+2 // /^+i {K^ - K^) ■ (ugV, (r+i - f+^) 

+ (/”+^ - V„ug) dvdx. 

(4.29) 

We will use the notation 

„n,l < 

to mean Ve > 03A^ such that n,i > N ^ cr”’* < e. Recall that 

sup||(T;-H-)-(if^i?0Li(M) + 

+sup// vS (/'*+! - f+^y dvdx < O"’^. 

By (4.14) we have 

\K'^-K^\ < Uo |(£;",H”) - 1^1 2-""/2 

^ l-j-e/2 n £ 

< Vq ' 

where we assume £ > n. 

Recall that a = a + 2 + e and note that 

// |/^+i (NT" - K^) ■ v^Vy (/"+! - f+^)\dvdx 

< |V„ (/"+i - f+^)jdvdx 

< //u“ dvdx \JII v^ |V^ (/”+^ - dvdx 

< d^’^yJJ v^ \Vy (/”+! - /^+i)|^ dvdx. 
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Similarly, 


// \f+^ {K^ - K^) ■ (/"+! - f+^) VyvS\dvdx 
< ff Vq (/"+^ — dvdx < o”’^. 


So, by (j4.29|) we find 

Tt II < (/"^^ - d^dx <-ff vS I V„ (r+i - f+^)l^dvdx 


+ 0 "’^ + 


0^’\/lIvS |V. (Z’^+i - /^+i)rdndx 


-//^o|V. (r+i-/^+i)|'dndx 

< o".^ - ff vg I (/"+! - dvdx. 

It follows that 

pT 


L I " f^^)\^dvdxdT < O^’^. (4.30) 


/o 

In a similar manner we may show that 

pT 


< |V„r+y dndxdr < C. 


(4.31) 


Note that the exponent of z;o in ()4.30|) is a, but in (j4.31ll it is a. 

Next, we derive bounds on Q. Considering 0 < t < t < T and letting 
p > 1, b,0 >0 and 


u = 


v — w X — y — K^{vi + wi) 

z = - 


y/t — T {t — t)3/2 

we have 

II WQ^^Q^dwdy 

= Cff (^1 + k-Vi-rwl") [(t - r)-5/2e-l“l V4e-3. 


{t - T)^/'^dzdu 

C{t - T)t(i-P) (^/|„|<+ e-Pl“l'/4du 

+ hi>i(t- , , , 


i|,|>i(t_,)~i/ 2 |„| (Vl + 1^^-Vi-rM|2) e 






(yrryp) 


-be 


+ e 


-C\v\^ 


< C{t - T)i^^-P'^Vo’‘'^. 
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So for 1 < p < 7/5 we have 


Wo^^gPdwdydT < . (4.32) 

Here, constants may depend on p, 6, and 9. Later, specific choices of p, b, and 9 are 
used and this dependence is removed. 

Next, we bound |Vu,C/| and |V„C/|. Note that 



If yj^bd 


= C 


II (^V^ + \v-y/t-Tu\y 


-be 


l“l 


'hit-r) 


',-5/2„-|«|V4„-3z^ 


= C{t - T)i-3P f (Vl +lv-^t- Tu\y 


-be r 


(t 


= -hlV4 


du 


- rf/^dzdu 


< C{t - 


and similarly. 



-y- + «^i)| 

{t-rY 


p 

dwdy < C{t 


t)2 


^-spy-be^ 


Hence, 


yy ZCO(iv^er + \yvg\P)dwdy < C(t- t)5-3Puo-'" 


and for 1 < p < 7/6 



tCo(|V„er + \V,g\P) dwdydT < 


In a very similar manner it may be shown that 



Vo^^gPdvdxdt < 


(4.33) 


(4.34) 


for p < 7 j5 and 



(|V„e|P + |V^e|P) dvdxdt < 


for p <11%. 

Next we derive two inequalities that will be used repeatedly. Let p,q € [l,oo) 
and r G [1, oo] satisfy 


1 

P 



We will first consider the case r Y oo, but what follows may be easily adapted to the 
case r = Qo. Let 0 > 0 and define 


b = 



c = 



-1 
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and note that 


111. 
- + - + - — 1 . 
roc 


dwdydr 
, 1/& 


l/c 


By Holder’s inequality and using (I4.32L we have for p < 7/5 and any h(T, y,w) >0 
fo If Qhdwdydr 

^ {fo If SP iwf)hy dwdydr^ ^ {f^ ff gPwg^^dwdyd? 

{fo If dwdydr'j ' 

< Cvo^ {fo ff {wohy dwdydr^ ' {/^ JJ (w^hy dwdyd-, 

Hence, using (|4.34ll we have, for p < 7/5 

fo ff (^0 fo ff Qhdwdydr^ dvdxdt 
fo ff {fo ff {w^hy dwdydrj 

{fo ff {'^o^y dwdydr^ dvdxdt 

{ly ff dwdydr^j ' 


< C 


< C 


Ijr 


< c 
= c 


1+r/c 

1/9 


.-I l/i- 


ly ff {fo ff dvdxdt^ {w^hy dwdydi 
{ly ff i^o^y dwdydr') 
ly ff {w^hy dwdydi 


l/r 


Similarly, using (I4.33[) we find for p < 7/6 


ly ff (^0 fo ff [l hdwdydr) dvdxdt 

^ ^ {ly ff {‘^o^y dwdydr') 


1 l/r 


(4.35) 


(4.36) 


Now we will show that /" converges in L°°. We have 

^„+i _ ^,+1 II . v„/^+i) (r, y, w)dwdydr. (4.37) 

Using (4.14) and taking i > n 

ITsT" . V.,/"+i -K^ -Vyf+^l < \K^-K^\ |V„/"+H 

+ \K^\ |v, {r+^-f+y\ 

< |(i?",H-)-(^;^H«)|^co|V„/"+l| (4.38) 

+C/u;o+^/"2-"^/2 I V„/"+i I + Cwo I V„ (/"+i - f+y \ 

< \v,r+^\ + cwo |v„ (r+i - f+y\. 
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Next, we will apply (14.351) with ip <1 jh and q = 2. Note that by taking p close to 
7/5, we may make r = (p + g~l) close to 


5 

7 



14/3. 


Thus, applying (14.351) with p < 1 jh, q = 2, 0 = | — 1, and 

h = «;o|V„(r+i-/+i)| 


while using (I4.30|) yields 


lo If (^0 Vo // ^^0 I V„ (/’"+! - f+^) I dwdydT^ dvdxdt 


- 1/r 


< C 


n{^l '«^o|v4/"+i-/'+i)|) dwdydT 


< o” 


for r < 14/3. Applying (14.351) with p < 7/5, q = 2, 9 = f — 1 — e/2, and h = 
|V„/"+^|, and then using (I4.31|) yields 


fo^ If ('^0 ^ dvdxdt 


1 1/r 


< c 


ji (»< 


f-l-e/2 l+e/2 


Wn 


v„/ 


n+1 


dwdydr 


< C 


for r < 14/3. Since f-l-f>f-l, (Him and (H3a now yield for r < 14/3 

rT 


/O 


f-1 

b 


r+^- 




dvdxdt < d" 


(4.39) 


Similarly, using (j4.31|) we may show that for r < 14/3 


— — 1 


ir-l)’ 


dvdxdt < C. 


(4.40) 


To use (14.391) we integrate by parts in (14.371) to obtain 


^n+1 _ _ 


• (if"/-+i - K^f+^) dwdydr (4.41) 


and using (4.14) 

\K^P+i - K^f+^ \ < - {E^,B‘)\wo |/”+^| 

I /"+! I + Cwo I /”+^ - I (4.42) 

< !/’"+!I + Cwo |/"+^ - . 

Next, we will apply (I4.36|) . but now with p < 7j6 and q < 14/3. Note that we 
may take p close to 7/6 and q close to 14/3 to make r close to 
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Thus, applying (14.361) with p <7 jd, q < 14/3, 0 = | — 2, and 


h = wo r+^-f 


n+1 


and then using (14.391) yields 


/o^ If (^0 ^ fo If \'^wQ\Wo II dwdydr^ dvdxdt 

1/9 


1 1/r 


< c 


j 7 is {< 


wo r+' - 


f+^lYdwdydi 


< d 


ni 


for r < 14. Applying (14.3611 again with p < 7/6, q < 14/3, 0 = ^ — 2 — | and 
/i = 

lo II lo II I dwdydrj dvdxdt 

< lo II \r^^\) dwdydrl <C 

for r < 14 by using (14.401) . With these two estimates, (14.411) and (j4.42|l yield 






dvdxdt < o" 


for r < 14. Similarly, 


ir-l)' 


dvdxdt < C. 


Finally, we can apply (14.361) with p = ■^,q = and r = oo. Proceeding as 
above we obtain 


,2 3 / rn+l fi+1 




< o”’F 


Recalling 


it follows that 


/= lim/" 

n—>-oo 


d dr-f) 


0 as n —?► oo. 


Next, we bound V„/" in L°°. We have 


fU + l 


GK^ ■ VyT+^dwdydT 


SO 


|v„ (r+i-iF)| = 
<cj;;n\^^g\wo 


lo II ^vGK^ ■ W^r+^dwdydT 

dwdydr. 


(4.43) 
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Applying (H3S1) with p < IjQ, q = 2, 6 » = f - 1, and h = wq |V„/"+i|, 
using (14.311) yields 


fo^ If ('^0 Vo // Wo I V„/"+^ I diudydr^ dvdxdt 

If |V„/”+^|) dwdyd 


- l/r 


1/2 


< c. 


Note that taking p close to 7/6 yields r close to 


J 2 

Hence, using (14.4311 we find for r < 14/5 


-1 


® + i - 1 ) = 14/5. 


/ 7 /(.r‘|v.r«l)’ 


dvdxdt < C. 


We then apply (I4.36p three more times. In each application we take h = wq 
First, using p < |, q < 14/5, and 9=^ — 2 yields 


• 7-2 




dvdxdt < C 


for r < 14/3. Using p < 7/6, q < 14/3, and 6 » = f - 3 yields 


f-3 




dvdxdt < C 


for r < 14. Using q = j. — 5 / = ^ _ 4 yields 


#-4 


Vo^ Vvf 


■n+1 


2 

< c. 


Recall that a> 8 so^ — 4 = — 4 > 1. Now V/i G 


/”+ {t,x,v + h) - /"+ it,x,v)\ < C\h\ 


and so 


\f{t,x,v + h) - f{t,x,v)\ < C\h\. 

Finally, we show that / is Holder continuous in x. Let h > 0 and 
e = P^^{t,x + h,v) - P+^{t,x,v), 

then 

dte + vidxG + iF" • V^e — A„e 

= — {K^{t, x + h,v) — K'^{t,x, v)) ■ X + h, v). 

Note that 

|Ar"(t, X + h,v) — K'^it, X, ?;)| 

< Vo (|A”(f, x + h)- E^{t, a;)| + x + h) - B'^{t, a:)|) 

< f {dxEf dx + f (dxBf dx^ 

< C'uohi/2. 


and then 


V„/"+i 


(4.44) 
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Thus, by (14.441) we find 


\dte + vid^e + ■ V^e - A„e| < C'uo/i^/^C'uo"“/^ < 

By Lemma 2.1 

|e| < 

It follows that / is Holder continuous with exponent 1/2 in x. Hence, Theorem H.l of 
[l6] applies and shows that / has the regularity stated in Theorem 1.1. The regularity 
of E and B follows from this. 

Finally, suppose that {F,£,B) is another solution with the same initial value as 
if,E,B). Then, by (glSD, dUS]), and dHH) 


Also 


Since 


^ II ~ II ~ 

+C J {\E - £f + {B - Bf) J E^v^o+^dvdx 

< C jj v^oif - Ffdvdx + C j {\E- £\^ + {B - Bf) dx. 


^ y (if; - + {B- Bf) dx = -2 j{E-£)- J if-F)v dvdx. 


(/ - I (f ~ J '^0 ^dv 

< cjif- F)^v^,dv 


(4.45) 


we have 


J {\E - £\^ - \B - B\^) dx < J \E- £\^ dx ^ JJ{f - F)\g dvdx 

< C J\E-£\^dx + C JJ{f-F)^v^odvdx. 

(4.46) 

Uniqueness follows from (14.451) and (14.461) and the proof is complete. 
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